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ABSTRACT 

We report on two new XMM-Newton observations of the Anomalous X-ray Pulsar 
(AXP) 4U 0142+614 performed in March and July 2004, collecting the most accurate 
spectrum for this source. Furthermore, we analyse two short archival observations 
performed in February 2002 and January 2003 (the latter already reported by Gohler 
et al. 2005) in order to study the long term behaviour of this AXP. 4U 0142+614 
appears to be relatively steady in flux between 2002 and 2004, and the phase-averaged 
spectrum does not show any significant variability between the four epochs. We derive 
the deepest upper limits to date on the presence of lines in 4U 0142+614 spectrum 
as a function of energy: equivalent width in the 1-3 keV energy range < 4eV and 
< 8eV for narrow and broad lines, respectively. A remarkable energy dependence in 
both the pulse profile and the pulsed fraction is detected, and consequently pulse-phase 
spectroscopy shows spectral variability as a function of phase. By making use of XMM- 
Newton and INTEGRAL data, we successfully model the 1-250 keV spectrum of 
4U 0142+614 with three models presented in Rea et al. (2007a), namely the canonical 
absorbed blackbody plus two power-laws, a resonant cyclotron scattering model plus 
one power-law and two log-parabolic functions. 

Key words: stars: pulsars: individual: 4U 0142+614 - stars: magnetic fields - stars: 
neutron - X-rays: stars 



1 INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) form a restricted family 
of neutron star (NS) X-ray sources with properties quite 
different from those of other known classes of X-ray emit- 
ting NSs, in particular the radio pulsars. However, in the 
last few years evidence gathered of a possible link between 
AXPs and another peculiar class of NSs: the Soft Gamma- 
ray Repeaters (SGRs; see Woods & Thompson 2006 for a 
recent review). 

The nature of the strong X-ray emission from AXPs 
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and SGRs is intriguing. In fact, their X-ray luminosity is too 
high to be powered by the loss of the star rotational energy 
alone, as in the more common radio pulsars. At the same 
time, the lack of observational signatures of a companion 
strongly argues against an accretion powered binary system, 
favouring instead scenarios involving isolated NSs. 

AXPs are usually radio-quiet and exhibit X-ray pulsa- 
tions with spin periods in the ~ 5-12 s range. This, along 
with their large spin-down rates {P ~ 10"^^ — lO^^^ss^^), 
implies magnetic fields exceeding the (electron) critical mag- 
netic field, Bcr ~ 4.4 X lO" G, if spin-down occurs mainly 
through magneto-dipolar braking. 

Furthermore, they are characterised by a high X-ray lu- 
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Figure 1. Top row: 0.3-10 kcV background-subtracted pulse profiles of all the four XMM-Newton observations of 4U 0142+614 . Bottom 
row: pulsed fraction as a function of energy for all the four XMM-Newton observations of 4U 0142+614 : the top and bottom panels 
report, respectively, on the pulsed fraction (as defined in the text) of the fundamental and first harmonic sine functions which best fit the 
pulse profile of each observation (note that first two and last two panels from left have different y-axis units). For both top and bottom 
rows, the relative observations arc (from left to right): February 2002, January 2003, March 2004 and July 2004. 



minosity {Lx ~ 10'''*-10''^ergs~^). Their high-energy emis- 
sion (~ 0.5-100 keV) is traditionaUy modelled by an ab- 
sorbed blackbody {kT ~ 0.4 keV) and two power-laws 
(rsoft ~ 3 and Fhard ^ 1), although recently several at- 
tempts have been made to apply more physical spectral 
models (e.g. Perna et al. 2001; Rea et al. 2007a, b, c). 

At present the model which appears most successful 
in explaining the peculiar observational properties of AXPs 
and SGRs is the "magnetar" model. In this scenario AXPs 
and SGRs are thought to be isolated NSs endowed with 
uhra-high magnetic fields {B ~ lO" - 10^^ G) the decay 
of which powers their X-ray emission (Duncan & Thomp- 
son 1992; Thompson & Duncan 1993, 1995, 1996). At mag- 
netic fields higher than Bcr, radio emission is believed to be 
suppressed e.g. by the photon-splitting process (Baring & 
Harding 1998), and this could explain why magnetars were 
believed not to show any radio emission. However, the dis- 
covery of radio pulsars with magnetic fields higher than the 
electron critical field (Camilo et al. 2000) and with very weak 
X-ray emission, as well as the discovery of radio pulsations 
from one magnetar (Camilo et al. 2006), is puzzling and still 
awaits explanation. Alternative scenarios, invoking accretion 
from a fossil disk, remnant of the supernova explosion (van 
Paradijs, Taam & van den Heuvel 1995; Chatterjee, Hern- 
quist & Narayan 2000; Perna et al. 2000; Alpar 2001), en- 
counter increasing difficulties in explaining the data. 

4U 0142+614 is one of the brightest AXPs known to 
date, and it was first detected by Uhuru in 1978 (Forman 
et al. 1978). However, mainly because of the presence of the 
accretion-powered binary pulsar RX J0146. 9+6121 nearby, 
only in 1994 was an ~ 8s periodicity reported using EX- 
OS AT data taken in 1984 (Israel, Mereghetti & Stella 1994). 



Long term spin period variations were discovered thanks to 
a large 7?XTE campaign (Gavriil & Kaspi 2002) , leading to 
the measure of the period derivative P ~ 2 x 10"^^ ss~^. De- 
spite deep searches (Israel, Mereghetti & Stella 1994, Wilson 
et al. 1999), no evidence for orbital motion has been found, 
supporting the isolated NS scenario. 

Subsequent observations (White et al. 1996; Israel et 
al. 1999; Paul et al. 2000) revealed a soft (1-10 keV) X- 
ray spectrum typical of an AXP, best fitted by an absorbed 
blackbody {kT ~ 0.4 keV) plus a power law (F ~ 3.7). 
More recent Chandra (Juett et al. 2002; Patel et al. 2003), 
XMM-Newton (Gohler et al. 2005) and Swift-XRT (Rea et 
al. 2007a) observations have shown that 4U 0142+614 is a 
relatively stable X-ray emitter. 

In the last few years, two peculiar characteristics of 
4U 0142+614, have been found in comparison with other 
AXPs: an optical counterpart (HuUeman, van Kerlwijk 
& Kulkarni 2000), displaying ~ 8s pulsation with a 30% 
pulsed fraction (Kern & Martin 2001 jj and li) mid-infrared 
emission, tentatively interpreted as the signature of a non- 
accreting disk around the NS (Wang, Chakrabarti & Ka- 
plan 2006). Furthermore, like in other AXPs, a hard X-ray 
emission up to 250 keV has been revealed (Kuiper et al. 2006; 
den Hartog et al. 2006, 2007). 

In this paper we report on two new XMM-Newton ob- 
servations of the AXP 4U 0142+614 performed in March 
and July 2004. These observations provide by far the most 



Note that the detection of the optical counterpart to this object 
might be due to its lower extinction with respect to other AXPs, 
rather than to an intrinsic peculiarity. 
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Figure 2. Epoch folding as a function of energy (from top to b( 
the 2004 March (left panel) and July (right panel) observations. 

accurate data collected in the soft X-ray band for this source. 
Furthermore, in order to study the long term X-ray activ- 
ity of 4U 0142+614 we analysed two archival observations 
(one of them already published by Gohler et al. 2005) . Data 
analysis and results are presented in § [5] while § |3] deals 
with the spectral modelling of the 1-250 keV spectrum of 
4U 0142-1-614 . Finally our results are discussed in §|4l 



2 OBSERVATIONS, DATA ANALYSIS AND 
RESULTS 

4U 0142-1-614 has been observed with XMM-Newton four 
times, the first two times for a very short exposure time 
(^6ks), while the last two observations are much deeper. 
These two very deep observations have been performed on 
March 1st and July 25th 2004, with on-source exposure 
times of 46 and 23 ks, respectively (see the end of this section 
for further details on the first two short observations). The 
XMM-Newton Observatory (Jansen et al. 2001) includes 
three 1500 cm^ X-ray telescopes with an EPIC instrument 
in each focus, a Refiecting Grating Spectrometer (RGS; den 
Herder et al. 2001) and an Optical Monitor (Mason et al. 
2001). Two of the EPIC imaging spectrometers use MOS 
CCDs (Turner et al. 2001) and one uses a pn CCD (Striider 
et al. 2001). 

Data have been processed using SAS version 7.0.0, and 
we have employed the most updated calibration files (CCF) 
available at the time the reduction was performed (Novem- 
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: 0.5-2 keV, 2-3 keV, 3-4 keV, 4-5 keV, 5-7 keV, 7-10 keV) for 

ber 2006). Standard data screening criteria are applied in 
the extraction of scientific products. We have cleaned the 
March 2004 observation for solar fiares (the observation in 
July is not affected); the exposure time after the solar flare 
filtering is 37 ks. 

Both 2004 observations have the same instrument set- 
up, pn and MOSl cameras are used with the meditmi filters, 
while the thick filter is applied for the MOS2 camera. The 
pn camera is set in timing mode in order to avoid pile-up, 
the MOSl with the central CCD in timing mode, and the 
others in prime full window mode. The M0S2 CCDs are 
all set in prime full window mode, with the aim of un- 
veiling possible transient X-ray sources which might con- 
taminate the instruments observing in timing mode. No 
transients are present in the MOS2 image. The source is, 
as expected, highly affected by pile-up. Hence, we do not 
consider the M0S2 data any further. 

In order to extract more than 90% of the source counts, 
we have constructed, for the data obtained with the pn and 
the MOSl central CCD, a one-dimensional image and we 
fit a Gaussian to the ID photon distribution. We have then 
extracted the source photons from a rectangular region with 
RAWX and RAWY coordinates 37,100.75 and 18,197.5. The 
background is obtained from a rectangular region of the 
same size, as far as possible from the source, on the east 
side. Only photons with PATTERN<= 4 are used for the 
pn and with PATTERN<= 12 for the MOSl. 

From the same regions we have extracted the 
4U 0142+614 spectra. All of the EPIC spectra have been re- 
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Table 1. Parameters of the spectral modelling of the 4U 0142+614 spectra with an absorbed blackbody plus a power-law, for all the 
four XMM-Newton observations. The last column reports on the parameters derived from a joint fit of all four spectra with the latter 
model. Uncertainties in this table are given at 90% confidence level. Nh is in units of 1022 cm~2 and solar abundances from Anders 
&; Grevesse (1989) are assumed. The blackbody radius is calculated assuming a distance of 3kpc (note that in the error calculation we 
did not consider the error in the distance). The blackbody fiux fraction is calculated with respect to the absorbed 0.5-10 keV total flux. 
Fluxes are all absorbed and in units of 10~io ergcm~2 . 



binned before fitting, so as to have at least 40 counts per bin 
and not to oversample the energy resolution by more than a 
factor three. We have also extracted first and second order 
RGSl and RGS2 spectra for both source and background, 
using the standard procedure reported in the XMM-Newton 
analysis manual. 

Thanks to the high timing and spectral resolutioijf] of 
the pn (0.03ms and E/(5E ^ 50) and MOS cameras (1.5ms 
and E/JE ~ 50), and to the high spectroscopical accuracy 
of the RGS (E/(5E ~ 400), we are able to perform timing 
and spectral analysis, as well as pulse-phase spectroscopy. 
The RGS spectra do not show evidence for narrow spectral 
features, and its continuum spectra are in good agreement 
with those of the pn and MOSl. We then do not report 
further on RGS data and refer to Durant & van Kerkwijk 
(2006a,b) for details on those. Furthermore, both MOSl and 
pn give consistent timing and spectral results. Hence, we re- 
port hereafter only on the pn data. We find a pn (background 
subtracted) count rate of 46.2 ±0.1 and 46.5 ±0.1 count/s 
for the March and July 2004 observation, respectively. 

Furthermore, we have also analysed two short (6 ks ex- 
posure time each) XMM-Newton observations performed 
on 2002 February 13th and 2003 January 24th (the lat- 
ter already reported in detail by Gohler et al. 2005). In 
both observations the pn camera is set-up in Small Win- 
dow mode. We use only pn data for these observations, and 
refer to Gohler et al. (2005) for details on MOS data for 
the January 2003 observation (MOS cameras gave results 
consistent with the pn for both observations). We have re- 
processed and cleaned the data as reported above for the 
2004 observations, resulting in an on-source exposure time 
in 2002 and 2003 of 2.0 and 3.5 ks, respectively. Thanks to 
the 2D image provided by the pn in Small Window mode, 
we have extracted the source events and spectra from a cir- 
cular region of 25" around the 4U0142±614 position. The 
background events and spectra, for both observations, are 
extracted from a similar circular region as for the source, 

2 see|http://xmm.esac.esa.int/|for details. 



although as far as possible (still within the central CCD) 
from the 4U0142±614 position. The source pn count rate 
is 42.1 ± 0.5 and 40.5 ± 0.6 count/s for the February 2002 
and January 2003 observations, respectively. 

2.1 Timing analysis 

In order to determine the spin period of 4U 0142+614 , we 
have first corrected the event times to the barycenter of the 
Solar System (using the SAS tool barycen) and computed 
a power spectrum (powspec), where we find a strong signal 
around 8.6 s in all the observations. We use for the timing 
analysis only events in the 0.3-10 keV energy range. Timing 
analysis is performed using Xronos. Two harmonics of the 
spin period are significantly detected in all four observations. 

A precise estimate of the pulse period in the observa- 
tions is obtained carrying out a period search (efsearch) 
around the 8.6 s value. We have then refined the period de- 
termination by dividing each observation into six intervals, 
and performing a linear fit to the phase determined in each 
interval (phase-fitting technique) . 

We obtain a best spin period of =8.6887(3) s (MJD 
52318.0), P.=8.6882(2)s (MJD 52663.0), Ps=8.68856(l) s 
(MJD 53065.0) and P,=8. 688575(2) s (MJD 53211.0), for 
the XMM-Newton observations from February 2002 to July 
20043. AH these periods are consistent within their 3cr errors 
with the more accurate RXTE timing analysis of this source 
(Dib et al. 2007). 

The pulse profile of 4U 0142+614 is double peaked in all 
observations, and well modelled by two sine functions with 
periods fixed at the spin period and its first harmonic (see 
Fig.[T} • Although the X-ray emission of 4U 0142+614 above 
100 keV is completely pulsed (pulsed fraction of ~ 100%; 
Kuiper et al. 2006), in the soft X-rays 4U 0142+614 shows 
the smallest pulsed fraction of the class (a factor of ~ 3 
below that observed in other AXPs). 

^ All errors in the text are reported at la confldence level, if not 
otherwise specifled. 
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The pulsed fraction (defined as the semi-amplitude of 
the best-fitting fundamental sine function divided by the, 
background corrected, constant level of the emission) in the 
0.3-10 keV energy band is 5.3 ±0.4%, 5.9 ±0.4%, 4.7 ±0.2% 
and 5.4 ± 0.2%, from February 2002 to July 2004. For a 
better comparison with previous results, we report also the 
peak-to-peak 0.3-10 keV pulsed fraction (defined as {Fmax — 
Frmn)/{Fmin + Fmax)), again fi^om 2002 to 2004: 8.8±0.9%, 
7.7 ± 0.9%, 10.2 ± 0.2% and 13.3 ± 0.3%. 

The 0.3-10 keV peak-to-peak pulsed fraction shows a 
modest (~ 2.5cr) increase between 2002 and 2004, possibly 
consistent with the evidence for a slow increase of the pulsed 
flux observed by 7JXTE (Dib et al. 2007). However, this is 
not followed by a similar increase of the fundamental sine- 
function pulsed fraction, which instead is within the errors 
relatively constant in time. This might be explained with 
the small changes of the pulse profile due mainly to the first 
harmonic of the signal, rather than to the fundamental. 

Studying in detail the pulsed fraction as a function of 
energy, in all the observations, we find a peculiar pulsed frac- 
tion increase with energy: e.g. in March and July 2004 the 
fundamental harmonic pulsed fraction increases from about 
4% to 15% as the energy increases from 1 to lOkeV (see 
Fig.[l} . Given the large errors on the pulsed fraction values in 
the 6-10 keV range for the February 2002 and January 2003 
observations, this increase of the pulsed fraction with energy 
is consistent with having the same trend (within the errors) 
in all four observations. This effect was already hinted in 
other observations but the lower number of counts collected 
was not sufficient to reliably asses it (see e.g. Fig. 5 of Israel 
et al. 1999, Paul et al. 2000, and Gohler et al. 2005). It is 
worth to note that this trend can not be an instrumental 
effect of the pn camera since the observations we report are 
taken in two different pn modes (Small Window and Tim- 
ing), and on the other hand this pulsed fraction increase 
was not observed in other AXPs or pulsars observed with 
the pn. Furthermore the presence of a hint of this increase 
also in BeppoSAX (Israel et al. 1999) and ASCA data (Paul 
et al. 2000) rules out any possible instrumental effect. 

Correspondingly, the pulse profiles are highly variable 
with energy, as shown in Fig.[2l This variability of the pulse 
profile with energy has always been observed in this source 
(Israel et al. 1999; Paul et al. 2000; Patel et al. 2003; Morii 
et al. 2005), as well as in other AXPs (see e.g. Rea et 
al. 2005; Mereghetti et al. 2004). To asses the statistical 
significance of these pulse profile variations, we have per- 
formed a Kolmogorov-Smirnov test on the energy resolved 
pulse profile for the March and July 2004 observations (see 
Fig-Ell 1 S'lid in both observations the difference between the 
soft (0.5-2 keV) and hard (7-10 keV) pulse profile is highly 
significant, giving a probability that the profiles have the 
same shape smaller than 10~®. We have done the same test 
between the two observations to asses the pulse profile vari- 
ability in time for a given energy range: in all the energy 
ranges the probability of the March and July 2004 profiles 
having the same shame is ~ 0.02, hence not very significant. 
However, note that RXTE monitoring over the last ten years 
did detect pulse profile variability with time, thanks to its 
higher sensitivity for pulse profile changes (Dib et al. 2007). 

This increase with energy of the pulsed fraction of the 
fundamental component has not been observed so far in 
other AXPs, while the pulse profile change with energy is a 
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Figure 4. Contours plot of the photon index (F) and the absorp- 
tion parameter for the joint fit reported in Tab. 1. Confidence 
levels are from small to large: la, 90% and 99%. The small cross 
labels the values of the joint fit parameters reported in Tab. 1. 

more common behaviour (see e.g. Rea et al. 2005; Mereghetti 
et al. 2004; Woods et al. 2004). 

2.2 Spectral analysis and pulse— phase 
spectroscopy 

The phase-averaged spectra for all the observations are sat- 
isfactorily fitted by an empirical model consisting of an ab- 
sorbed blackbody plus a power law (Xdo/ — 1-0; see Tab. 1 
and Fig. [3]). We use the XSPEC package (version 11.3 and 
12.1) for all fittings, and used the phabs absorption model 
with the Anders & Grevesse (1989) solar abundances. A 
2% systematic error is included in all fits, as well as in 
the error calculation, to account for calibration uncertain- 
ties. Similarly to what has been done in the past for other 
AXPs, we have tried to fit the spectra of 4U 0142+614 with 
two absorbed blackbodies and with a Comptonized black- 
body, but both trials resulted in a worse fit (Xdo/ > 2.4). 
In all the observations the hydrogen column density is 
Nh ~ 1 X lO^^cm"^, and the absorbed flux in the 0.5- 
lOkeV range is ~ 1.1 x 10^^" ergcm^'^ s~^, corresponding 
to a luminosity of ~ 1 x 10^^ ergs~^ (for a 3 kpc distance). 
In the 0.5-10 keV band, the blackbody component accounts 
for ~ 24 % of the total absorbed flux. The blackbody radius, 
as derived from the fits, is smaller than the NS size. If the 
blackbody emission originates from the star surface, and is 
then partially up-scattered to form the power-law compo- 
nent as predicted by the twisted magnetosphere model (see 
33] and Q , this would imply that only a fraction of the sur- 
face is emitting. Q. 

^ Note that the measure of the emitting region can be strongly 
affected by the choice of the spectral model. For instance, in the 
case of cooling neutron stars, accounting for atmospheric effects 
would result in an emitting area larger than that inferred assum- 
ing blackbody emission (Perna et al. 2001). These atmospheric 
models are not direct applicable to AXPs, which are unlikely to 
be passive coolers. Nevertheless, the uncertainty on the spectral 
shape for the thermal component affect the determination of the 
emitting size. 
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Figure 3. Phase average spectra of the March (left) and July (right) XMM-Newton observations of 4U 0142+614 fitted with an 
absorbed blackbody plus a power— law model. 



The best fitting parameters for the absorbed blackbody 
plus power law model are only marginally (< 2a) variable 
between the observations (see Tab. 1; especially between the 
March and July 2004 pointings), and are all consistent with 
previous measurements (Israel et al. 1999; Paul et al. 2000; 
Juett et al. 2002; Patel et al. 2003). In order to shed hght 
on possible marginal spectral variability in the 2004 obser- 
vations, we have fit all the spectra simultaneously with an 
absorbed blackbody plus power-law model, with all parame- 
ters tied between the spectra. Results are reported in Tab. 1 
last column, and show that the spectrum does not display 
significant variability between the four observations. We also 
report in Fig.|3]on the contours plot for the photon index and 
absorption parameter for this the joint fit, showing that both 
parameters are relatively well constrained by this modelling. 
However, looking at the residuals we see a (not significant) 
excess in the residuals above 8keV (see e.g. FigO, proba- 
bly due to the onset of the hard X-ray component (Kuiper 
et al. 2006; den Hartog et al. 2006, 2007). To shed light on 
this issue, we have attempted the fit of the March 2004 ob- 
servation with the addition of the 25-250 keV INTEGRAL 
spectrum (see ij3]for further details). 

We also find some (not significant) deviations from the 
best fit model in the lower energy range (see Fig.|3]) for 
the two deepest observations in 2004. However, we believe 
these deviations are due to both remaining calibration issues 
and/or, possibly, deviations in the ISM abundances from 
the (assumed) solar ones (see also Durant & van Kerkwijk 
2006a). 

We have performed a pulse-phase spectroscopy analy- 
sis for all the observations, but we report below only the 
two deeper observations in 2004, because previous XMM- 
Newton observations gave consistent results but with a 
smaller accuracy (see also Gohler et al. (2005) for the PPS 
of the January 2003 observation) . 

We have generated 10 phase-resolved spectra for each 
of the 2004 observations. Phase zero is arbitrarily set close 
to the start of the observations (MJD = 53065.38674 and 
MJD = 53210.30000 for the March and July observations, 
respectively) and with all phase bins of the same size. The 
choice of the number of intervals has been made a priori 
in order to have enough statistics in each phase-resolved 
spectra to detect, at a 3a confidence level, a cyclotron line 



with an equivalent width > 30 eV. This allow us to keep 
the number of trials to a minimum in case evidence for a 
cyclotron line is found. 

The absorbed blackbody plus power law model provides 
an excellent fits for all the ten pulse phase-resolved spectra 
in both observations, either when leaving Nh free (see the 
2nd, 3rd and 4th panels from top in Fig.[S| or when fix- 
ing it to the value obtained from the phase-average anal- 
ysis (see last two panels in Fig.[5ll. In all the observations 
there is a similar spectral variability with phase, hence we 
show in Fig. [5] only the longer March 2004 observation. The 
phase-dependent photon index (as obtained from the fits 
with fixed absorption, see last panel of Fig.[SJ shows devi- 
ations from a constant value with a significance of about 
4(7, while the blackbody temperature shows a 3.3a variabil- 
ity. The blackbody temperature seems to vary with phase in 
anti-correlation with the photon-index, as also suggested by 
Patel et al. (2003). Note that Gohler et al. (2005) did not 
detect any phase-resolved blackbody temperature variability 
because of the much lower number of counts. 

However, it should be noticed that the Nh value we 
obtain from this modelling might be overestimated, since a 
lower value was found from fitting the ISM edges directly in 
the RGS data (Durant & van Kerkwijk 2006a) and by using 
other more physical, although yet less common, models (see 
33 and Rea et al. 2007a). 

No spectral features are detected in the phase-averaged 
and phase-resolved spectra of both the pn and RGS data 
(see also Juett et al. 2002 for the Chandra grating upper 
limits) of the two 2004 deep observations. In Tab. 2 we re- 
port the 3a upper limits of the equivalent width of a Gaus- 
sian line with o-iinc=0 (narrower than the pn instrumental 
energy resolution) and CTunc^lOOeV. These are the deepest 
upper limits to date on the presence of lines in the spec- 
tra of an AXP. So far, the only evidence for an absorp- 
tion line in an AXP, detected at a significance of ~ 3a, 
was discovered by BeppoSAX (Rea et al. 2003; but see 
also Rea et al. 2005, 2007d) in an observation of the AXP 
IRXS J170849.0-400910 . This absorption line, tentatively 
interpreted as due to resonant cyclotron scattering, was ob- 
served at an energy of 8.1 keV, with a width of 0.2 keV and 
an equivalent width (EW) of 0.8 keV. For IRXS J 1 70849.0- 
400910 , the absorption line was observed when the source 
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Figure 5. Top panel reports on the 0.3-10 keV pulse profile for 
the March 2004 observation. From top to bottom: 2nd, 3rd and 
4th panels report on the spectral parameters of the pulse phase 
spectroscopy analysis for the 4U 0142+614 March 2004 observa- 
tions (see also §[2]2j with the absorption parameter Nj{ free to 
vary. Last two panels correspond to phase-resolved spectral pa- 
rameters with fixed A^^f = 1 X lO^'^cm"^. Vertical dashed line is 
at phase 0.5 for clarity. 

had a peculiar high X-ray flux and hard spectrum, while 
was not detectable during a subsequent observation (two 
years later) when IRXS J170849.0-400910 was found at a 
lower flux level and with a softer X-ray spectrum (Rea et 
al. 2005, 2007d). Comparing our upper limits on the presence 
of lines in the spectrum of 4U 0142-1-614 (Tab. 2) with the 
IRXS J170849.0-400910 line properties (Rea et al. 2003), 
we could have easily detected in the 4U 0142-1-614 0.3-9 keV 
spectrum any absorption or emission line with the same 
width and EW as observed in IRXS J170849. 0-400910 . 



3 MULTIBAND SPECTRAL MODELLING 

The 4U 0142+614 X-ray spectrum has been discovered to 
have a hard X-ray component (Kuiper et al. 2006; den Har- 
tog et al. 2007). Rea et al. (2007a) applied several models to 
the 0.5-250 keV 4U 0142+614 spectrum in order to model 
the emission of this source across this large energy range, 
using simultaneous observations performed in July 2005 
with Smft-XRT (0.5-8 keV) and INTEGRAL (25-250 keV; 



Energy Range 


o-linc = (eV) 


o-iinc = 100 (eV) 


1-2 keV 


< 4 


< 9 


2-3 keV 


< 4 


< 7 


3-4 keV 


< 12 


< 19 


4-5 keV 


< 16 


<21 


5-6 keV 


< 25 


< 37 


6-7 keV 


< 46 


< 51 


7-8 keV 


< 88 


< 124 


8-9 keV 


< 405 


< 416 


9-10 keV 


< 527 


< 789 



Table 2. 3fT upper limits on the equivalent width of a Gaus- 
sian line with width o-iinc=0 (narrower than the pn instrumen- 
tal energy resolution) and 100 eV, derived from the March 2004 
XMM-Newton spectrum of 4U 0142+614 . 

den Hartog et al. 2007). Here we use the latter INTEGRAL 
spectrum fitting it together with our March 2004 XMM- 
Newton observation (being the longest one). This multiband 
fitting is performed in order to i) asses the presence of an 
hard component in the 8-10 keV XMM-Newton spectrum of 
4U 0142+614 (see Fig[3}, and ii) to test whether the spec- 
tral parameters derived in Rea et al. (2007a), are consistent 
with those derived from the higher XMM-Newton statistics 
and larger energy band (we used here the pn in the 1-10 keV 
energy range). 

We note that, at variance with what has been done in 
Rea et al. (2007a), these observations are not simultaneous. 
However, the soft X-ray spectrum shows only a very modest 
variability between the Swift-XRT ohserv&tion performed in 
July 2005 (simultaneous to the INTEGRAL observation we 
use here) and the XMM-Newton observations we report 
here. We then tentatively rely on the stability of the hard 
X-ray spectrum and try to fit the overall spectrum assuming 
that the source is in the same emission state during both the 
XMM-Newton and the INTEGRAL observations. 

We fit the 1-250 keV spectrum with: one absorbed 
blackbody plus two power-laws, an absorbed resonant cy- 
clotron scattering model plus a power-law and two absorbed 
log-parabolae (Tab. 3; see Rea et al. 2007a for further details 
on these models). 

We find that all these models, remove the excess in 
the >8keV residuals (not significant though) found in the 
XMM-Newton spectrum when fitted with an absorbed 
blackbody plus a single power-law (see Fig [3] and Tab. 1). 
We then conclude that this excess in the residuals is indeed 
due to the presence of the hard X-ray component. Further- 
more, the spectral parameters we find for the three best fit 
models are consistent (within Sct statistical accuracy) with 
those reported by Rea et al. (2007a) using the Swift-XRT 
spectrum, with lower statistics and smaller energy coverage. 
We note that the constant parameter in the fitting of the 
two log-parabolae is slightly smaller than the range of val- 
ues that the cross-calibration between XMM-Newton and 
INTEGRAL should take. This might be either a hint of 
the fact that the second log-parabola can barely take into 
account the 8-10 keV part of the spectrum which was miss- 
ing in the Swift-XRT spectrum, or simply a statistical fluc- 
tuation. A simultaneous observation of 4U 0142+614 with 
XMM-Newton and a longer INTEGRAL exposure might 
shed hght on this issue. 
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BB + 2PL 



RCS + PL 



2 log-pajrabolae 



Nh 

constant 
kT (kcV) 
BB Flux 

Tsoft 
^hard 

PLsoft Flux 

PLhard Flux 



926+° °°^ 

1.2 

422+° °°^ 
3.86«:«1 

'•"-1.0 
I 7+0.4 
^•'-1.2 



kT (keV) 
To 

13th 

Thard 

RCS Flux 
PLhard Flux 



0.51 



0.334 
1.84 
0.22 
1.1 



+0.01 
■0.01 

0.86 

+0.001 
-0.002 
+0.07 
-0.03 
+0.01 
0.02 
fO.l 
-0.1 
9 ,+0.3 
^•^-0.7 
1 r;+0.8 
^•^-0.9 



O.SS"* 



Epi (keV) 

01 

logPi Flux 
Ep2 (keV) 

ZosP2 Flux 



rj+O.Ol 

-0.01 

0.6 

, 07+0.03 

^■•^'-0.04 
+0.12 
-0.15 
+1.0 
-0.9 

200^8° 

-0 6+°-l 
"■"-0.1 

2 0+1-2 



-3.21 
2.1 



Total Abs. Flux 

Total Flux 
Xl (d.o.f.) 



^■°-0.5 

8 4+0-9 
^■^-0.3 

1.09 (225) 



3.0 



+0.1 



^-0.4 

1+0.8 
-0.5 

1.06 (225) 



4.r 



5+0.2 

-1.0 
hi.i 
-1.0 
0.82 (225) 



3.3: 

4.4"* 



Table 3. Best fit parameters for the 1—250 keV spectral modelling of 4U 0142+614 . Errors are at 90% confidence level. Fluxes (if not 
otherwise specified) are unabsorbed and in units of 10"^" ergcm"^ s~^. RCS, BB, PLgoft and logPi fluxes are for the 0.5-10 kcV energy 
band, while PLhard and logPi fluxes refer to the 20-250 keV band. Total fluxes are in the 0.5-250 keV band. Nu is in units of lO^^-^ cm"^. 
The constant parameter, which accounts for the inter-calibration, assumes XMM-Newton as a reference. 



4 DISCUSSION 

Wo reported here on the timing and spectral analysis of 
four XMM-Newton observations of the AXP 4U 0142-1-614 
performed from 2002 to 2004, representing the deepest X- 
ray observations on this source performed to date. 

The pulse profile of 4U 0142-1-614 is double peaked in 
the soft X-rays, and largely variable as a function of energy; 
we also detect an increase with energy of the fundamental 
pulsed fraction component (see also hints for this increase 
in Israel ct al. 1999; Paul ct al. 2000; Gohlor ct al. 2005). 
Furthermore, ten years of RXTE monitoring also revealed 
the pulse profile to be very variable in time, possibly related 
to the presence of glitches (Dib et al. 2007; but see also Morii 
et al. 2005). 

The presence of two distinct peaks in the phase profile 
of 4U 0142-1-614 suggests that the emission might be con- 
centrated in two regions on the star surface. Since the two 
maxima are offset by half a phase cycle, the emission re- 
gions must be located antipodally on the surface of the star. 
In a "standard" cooling NS, such a geometry of the emis- 
sion regions would be naturally produced by a (core-centred) 
dipolar magnetic field. In fact, since heat flows preferentially 
along the magnetic field lines, in this scenario the polar re- 
gions (where the field lines are more closely packed) turn out 
to be hotter than the equatorial belt. On the other hand, 
for such configurations, the pulsed fractions are low, the 
more when gravity eflfects are accounted for (e.g. Page 1995; 
Perna, Heyl & Hornquist 2000; DeDeo, Psaltis & Narayan 
2001). The very high pulsed fraction of the hard (> 100 keV; 
Kuiper et al. 2006; den Hartog et al. 2007), photons may be 
explained if they are produced in regions high up in the mag- 
netosphere, where no relativistic effects are expected to wash 
out the pulsed emission. This might also tentatively explain 
the increase of the pulsed fraction with energy, in particu- 
lar between the band dominated by the thermal component, 
and that dominated by the power-law. However, much more 
detailed theoretical models are need to fully interpret this 
phenomenon. In fact, it is unlikely that magnctars are pas- 
sive coolers in which the surface temperature distribution 
is solely determined by the large-scale field topology. As 
stressed by Thompson, Lyutikov & Kulkarni (2002), por- 



tions of the star surface are heated by the returning currents 
which flow into a "twisted" magnetosphere, and this may be 
a possible explanation for the presence of hotter regions on 
the star surface. It has been recently suggested (Thompson, 
Lyutikov & Kulkarni 2002; Lyutikov & Gavriil 2005; Fer- 
nandez & Thompson 2007) that the mechanism responsible 
for the formation of the 0.1-10 keV spectrum of AXPs and 
SGRs is cyclotron resonant up-scattering of soft (~ 0.5 keV) 
thermal photons produced at the star surface by electrons 
(or pairs) which fill the (twisted) magnetosphere. Currents 
are mainly concentrated towards the magnetic equator, and 
this, as the Monte Carlo simulations of Fernandez & Thomp- 
son (2007) show, makes the emission pulsed even if primary 
thermal photons are produced uniformly at the star surface. 
We performed some preliminary calculations with a Monte 
Carlo code (Nobili et al. in preparation) in order to check 
if the increase of the pulsed fraction with energy is quali- 
tatively reproduced by the resonant scattering model. For 
instance, assuming that the star is an orthogonal rotator 
seen at an angle ^ 90° with respect to the rotation axis, 
we find that the low energy (0.5-2 keV) light-curve should 
be double peaked while, at higher energies (7-10 keV), the 
secondary peak tends to disappear and the pulsed fraction 
increases. 

At variance with what has been recently found for a few 
other AXPs (e.g. IRXS J170849. 0-400910 , IE 1048.1 5937 
and IE 2259+586 ; Rea et al. 2005; Mereghetti et al. 2004, 
and Woods et al. 2004), these very deep X-ray observations 
of the AXP 4U 0142-1-614 find the source stable in flux 
and spectral shape. However, spectral and flux variability 
arc expected to be connected with the AXP bursting ac- 
tivity, as e.g. detected for IE 2259-1-586 (Kaspi et al. 2003; 
Woods et al. 2004). Hence, the recent bursting activity of 
4U 0142-1-614 (Kaspi et al. 2007; Gavriil et al. 2007) should 
have probably disturbed its steady emission state, and the 
source might have recently enhanced its flux and hardened 
its spectral shape. 

Furthermore, we put very strong upper limits on the 
presence of X-ray lines in the spectrum of this AXP, either of 
resonant cyclotron or atmospheric nature, the deepest upper 
limits for an AXP to date. This is a further confirmation that 
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the X-ray spectral lines, possible due to proton cyclotron 
scattering, or to absorption by the NS atmosphere, are either 
very weak, or not present at all (e.g. suppressed by QCD 
processes; Ho & Lai 2003; van Adelsberg & Lai 2006) in the 
spectra of AXPs, and probably also of SGRs. 



This paper is based on observations obtained with 
XMM-Newton , an ESA science mission with instruments 
and contributions directly funded by ESA Member States 
and the USA (NASA). NR is supported by an NWO Post- 
doctoral Fellowship, and thanks Lucien Kuiper for useful 
comments, and Vicky Kaspi and Marjorie Gonzalez for use- 
ful discussion about the pulsed fraction variability of this 
source. SZ thanks PPARC for financial support. We thank 
the referee for his/her very useful suggestions. 
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